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Long-lived single-frequency oscillations in the local non-equilibrium dynamics of a quantum many-body sys¬ 
tem is an exceptional phenomenon. In fact, till now, it has never been observed, nor predicted, for the physically 
relevant case where a system is prepared to be quenched from one quantum phase to another. Here we show how 
the quench dynamics of the entanglement spectrum may reveal the emergence of such oscillations in a corre¬ 
lated quantum system with Kondo impurities. The oscillations we find are characterized by a single frequency. 

This frequency is independent of the amount of energy released by the local quench, and scales with the inverse 
system size. Importantly, the quench-independent frequency manifests itself also in local observables, such as 
the spin-spin correlation function of the impurities. 


I. INTRODUCTION 


The dynamics of an isolated quantum many-body system 
following an instantaneous change of a Hamiltonian param¬ 
eter (quantum quench) is a topic of growing interest^. The 
problem touches on a multitude of subjects, from the founda¬ 
tion of quantum statistical physics to the engineering of quan¬ 
tum states and devices^il. A quench injects energy which dis¬ 
perses among the interacting degrees of freedom, and as time 
evolves local observables relax to their equilibrium values^. 
When the quench is local, with a sudden change of a local 
parameter in a Hamiltonian, the energy injected to the system 
is nonextensive. In such a case one may expect, and indeed 
finds^, intermediate times at which wave propagation and re¬ 
flection from boundaries can create slowly decaying oscilla¬ 
tory behaviour. This raises the question whether one could 
And a physically relevant model where, by tuning a pertinent 
parameter, the equilibration after a local quench of the Hamil¬ 
tonian is strongly suppressed, or even totally eliminated. 


Here, in point of fact, we show that, for the spin emulator— 
of the two-impurity Kondo modeli^, a local quench into the 
Kondo-screened phased induces the onset of long-lived os¬ 
cillations. By analyzing the quench dynamics of the lowest 
eigenvalues of the entanglement spectrumi^dS, we And that 
the frequency /„ of these oscillations is sharply determined 
if one tunes the impurity-spin interaction so that all spins be¬ 
come entangled with an impurity. Remarkably, /„ is indepen¬ 
dent of the local quench energy and scales as fuN remaining 
constant when the system size N increases. Moreover, we find 
that the frequency /„ leaves its fingerprints on local observ¬ 
ables, including the spin-spin impurity correlation functions. 
As such, we believe that our finding is potentially relevant to 
future designs of spin-based quantum devices as the impurity 
spins do not seem to equilibrate despite being strongly con¬ 
nected to reservoirs. 
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FIG. I: Schematic picture of the spin emulator of the two-impurity 
Kondo model. The impurity coupling J', together with the RKKY 
coupling K, determines the extension of the entanglement length 
^{J',K). The strength of the RKKY coupling K can be externally 
controlled; for instance, in a quantum dot, by tuning gate voltages 
as indicated in the figure. For any K < Kc, in the Kondo regime, 
we define as the impurity coupling for which the entanglement 
length extends over the entire chain, i.e. ^{Jopt,K) ~ N. By 
quenching the RKKY coupling K from Ki in the RKKY regime 
to K 2 in the Kondo phase one induces the non-equilibrium dynamics 
being studied here. For the sake of simplicity the next-nearest cou¬ 
pling J 2 is not shown in the figure since we fixed J 2 /Ji = 0.2412. 


II. MODEL 

We consider a spin chain emulator of the two-impurity 
Kondo model^ with two localized spin-1/2, each coupled to 
a frustrated spin-1/2 Heisenberg chain, and to each other via 
a Ruderman-KitteTKasuya-Yosida (RKKY) interaction. The 
Hamiltonian can be written as 


H= Y. {HZp + HZk} + Hi, (1) 

m—L,R 
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i—2 i—2 

Hi = KJicr^-o-^. (2) 

Here m = L,R labels the left and right chains, with cr™ the 
vector of Pauli matrices at site i on chain m. The couplings 
Ji and J 2 are nearest and next nearest neighbour couplings, 
respectively. The dimensionless parameter J' > 0 plays the 
role of an antiferromagnetic Kondo coupling between the im¬ 
purities and their corresponding bulks and the dimensionless 
coupling K > 0 measures the RKKY impurity interaction. 
The mapping to the Kondo model is valid for J 2 < JI (with 
J| = 0.2412 since in this interval the bulk excitations 

are massless (as in the full electronic version of the Kondo 
problem). As we will see, our analysis of the quench dynam¬ 
ics applies to the entire Kondo regime, although, for J 2 < J| a 
marginal coupling (in the sense of the renormalization group) 
produces logarithmic corrections which pollutes the numer¬ 
ical dat a*^'*^ . In order to avoid this, we tune J 2 = J| for 
which both Hl and Hu faithfully represent the spin sec¬ 
tor of a single-impurity Kondo model. Note that the terms 
ill E*!- O define a spin emulator of the 
single channel Kondo model, with labels m = L, A 
schematic picture of the model is given in Fig.[T] 

The ground state shows a quantum phase transition at a crit¬ 
ical value K = Kc of the RKKY coupling. For a small cou¬ 
pling K the Kondo interaction dominates and each impurity 
spin gets screened by its bulk (Kondo phase), while in the op¬ 
posite limit the impurities form a local singlet (RKKY phase) 
and decouple from the rest of the system^. 

III. RESULTS 
A. Entanglement spectrum 

Having partitioned the system into two parts, L and R (see 
Fig.ffl), one can write an arbitrary pure state in the orthogonal 
Schmidt basisi^ as 

\fpLR) = y/^lLn) 0 \Rn), > 0, (3) 

n 

where the ordered set of real numbers Ai > A 2 > ... form the 
entanglement spectrum. The Schmidt bases {| and 

{\Rn)}nJi diagonalize the reduced density matrices pL and 
pn of the left and right parts respectively. It has been shown 
that the two largest Schmidt numbers dominate the entangle¬ 
ment spectrum, with their difference, Ai — A 2 , behaving as an 
order parameter at the quantum phase transitio n* A Notably, 
all levels of the entanglement spectrum contribute essentially 
to the von Neumann entropy 

s{pl) = s{pr) = - ^ A„ log(A„). (4) 

n 

In the Kondo regime K < Kc, the system supports a length 
scale J', K) which diverges at the critical point in the ther¬ 
modynamic limit. It may be determined numerically by ex¬ 
ploiting its interpretation as the length scale over which the 


two impurities are entangled with two identical blocks of spins 
on both sides^. One thus finds, for a large but finite system 
with N = Nr 

^ g-a/j' 

iiJ \K - Kc\‘' + Oil/N)' 

Here a is a constant, with r a critical exponent taking the 
value v = 2 in the neighborhood of the critical point Kc ~ 
. According to Eq. (|5]), for any given K < Kc in 
the Kondo regime, it is always possible to find an optimal im¬ 
purity coupling J' = J'cpi such that ^{Jcp^,K) = N, mak¬ 
ing each impurity entangled with all spins in its bulk. For 
J' <C J'opt the length ^ exceeds the length of the chain and 
Kondo screening does not take placed. 

The time evolution of entanglement spectra following a 
quantum quench have been the subject of several recent 
investigations^^— . As we shall see next, the entanglement 
dynamics of the present problem exhibits some striking fea¬ 
tures. 


B. Entanglement dynamics 

We initially prepare the system in the ground state 
|GS(Ari)) of the Hamiltonian in Eq. ([T]i, choosing the cou¬ 
pling K = Ki > Kc in the RKKY regime. At time t = 0 the 
coupling is instantaneously changed to K = K 2 < Kc- As a 
result, the system evolves as 

n 

where |4'(0)) = |GS(Ari)), and where {i?„} are eigenvalues 
of the quench Hamiltonian, defined by Eq. (1T]| with K = K 2 . 
The corresponding eigenstates {|5'n)} are global singlets, as 
implied by spin-rotational symmetry. By tracing part R (see 
Fig.D, one can compute the reduced density matrix pL{t) of 
part L, from which the entanglement spectrum is obtained. 

Eigs. W[a) and (b) show the results for two quantum 
quenches with different values of Ki using exact diagonaliza- 
tion, with the two largest Schmidt numbers Ai (t) and A 2 (f) 
plotted as functions of time. In both cases J' = J'^pt, so that 
^{J'opt, K 2 ) = N. While the amplitudes of the oscillations are 
different for the two quenches, very surprisingly, the dynam¬ 
ics of A 1.2 is governed by a single frequency /, independent 
of the quench. Moreover, as seen in Eigs. |2|a) and (b), the 
oscillations do not damp. 

It is instructive to study the power spectrum 

A„(/) = J-[A„(f)] = ^ J n= 1,2. (7) 

In Eig. EJc), Ai(/) is plotted for the two given quenches. 
As revealed by the plot, a sharp peak, independent of the 
quench, emerges at the frequency f = fu- This shows that 
there exists a a unique frequency fu excited by quenching. In 
Figs.|2jd)-(f) the same quantities are plotted for a larger chain, 
using a time-dependent density-matrix renormalization group 
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FIG. 2: Quench dynamics of the first two largest Schmidt numbers Ai and A 2 versus time for J' = J'^pt = 0.3 in a chain of length N = 20 
with K 2 ~ 0.19, and Ki = 1 (blue solid line); K\ = 0.25 (red dashed line) using exact diagonalization: plots (a) and (b). The Fourier 
transform Ai (/) for J' = J'^pt = 0.3 in a chain of length N — 20 with K 2 = 0.19 and Ki = 1 (blue line with small solid circles representing 
data points) and Ki = 0.25 (red line with stars representing data points): plot (c). Quench dynamics of the first two largest Schmidt numbers 
Ai and A 2 versus time for J' = J'^pt = 0.24 in a chain of length = 40 when Ki = 1 and K 2 = 0.08 using tDMRG: plots (d) and (e). The 
Fourier transform Ai (/) for J' = J'„pt = 0.24 in a chain of length = 40 with Ki = 1 and K 2 = 0.08: plot (f). Quench dynamics of the 
first two largest Schmidt numbers Ai and A 2 versus time for a non-optimal case, J' 7 ^ (here J' = 0.7), in a chain of length N = 20 when 
K\ — 1 and K 2 = 0.19: plots (g) and (h). The Fourier transform Ai (/) for J' = 0.7 7 ^ J'opt in a chain of length A = 20 with K\ = 1 and 
K2 = 0.19: plot (i). 


(tDMRG) algorithm^^. Again, single-frequency oscillations 
are uncovered on the numerically accessible time scale. To a 
high precision, can be fit by a sinusoidal function as 

A„ = A„sin(27r/„f+ ^„), n = l,2 (8) 

where fu = {E 2 — Ei)/2'!t is the single frequency involved 
in the dynamics, with i?i 2 the eigenvalues of the states | 5 'i^ 2 ) 
which dominate the expansion in Eq. (| 6 ]l. and are the 
quench-dependent amplitude and phase. 

Of course, the single-frequency oscillations disappear for 
non-optimal Kondo coupling J' 7 ^ as shown in Figs. |2lg) 

and (h). Indeed, these plots clearly display dispersive multi¬ 
frequency oscillations, also evidenced by the power spectrum 
in Fig. |2i). It is, however, very important to observe that 
there always exists a dominant frequency fpeak in the power 
spectrum of Ai 2 which exhibit a pertinent scaling with J'. It 
follows from Fig.^a) that fpeak scales with J' as the Kondo 


temperature, 

fpeak ^ Tk, (9) 

using that T/f ~ 1/^{J',K2) ~ exp(—a/ J')i. For J' = 
J'p(, i.e. when ^(J'p(,A 2 ) = N, this frequency becomes 
the only one accessible to the system, i.e. fpeak = fu- For 
J' < J'pj we see deviations from Eq. (|9]l due to finite-size 
effects, as now J', K 2 ) exceeds the length N. 

It is important to realize that the single-frequency oscil¬ 
lations unveiled by our simulations are very different from 
the quench dynamics studied in Refs 4 i 2 ii 2 i, which reflects a 
finite-size effect when instantaneously joining two quantum 
critical systems. Likewise, the oscillating quench dynamics 
numerically observed in certain ID lattice models^Sril are also 
different from our results, coming from a global quench of an 
integrable interaction, and being either damped^^ or exhibit¬ 
ing a multi-frequency power spectrum^Sill. In all these cases 
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one observes dispersive wave propagation while our scenario 
is nonperturbative, long-lived and dispersionless. 


C. Quench independence 

The frequency /„ is quench-independent as it depends nei¬ 
ther on Ki nor 7^2- Independence from is evident from 
Eq. Independence from K 2 comes about by tuning the 
Kondo coupling J' to its “optimal” value J'pt, thus entan¬ 
gling the impurities with all spins in their respective bulks, 
i.e. A" 2 ) = N. For instance, in a chain of length 

N = 20, for K 2 = 0.19 and = 0.30 one finds that 
fu = 0.11; the very same value of /„ is obtained for K 2 = 0.1 
and J'p( = 0.315. 

The dependence of fu on the optimal Kondo coupling 
is plotted in Fig. Ob), showing that 

( 10 ) 

as expected from Eq. (O. By combining Eq. (fTOl) with Eqs. 
(|5]l and (O, it follows that for a chain of arbitrary size N, 

ajopt,K2) ^N' 

To check Eq. (fTTT i. we have computed the product fuN for 
various lengths N, keeping K 2 fixed; see TABLE 1 which 
shows that fuN ~ 2.1 ± 0.15. The scaling in Eq. (fTTIi sug¬ 
gests that we get data collapse onto a universal curve for dif¬ 
ferent lengths N if plotting Ai 2 vs. fut. This is confirmed in 
Fig.Oc). 


D. Effective model for J' = . - 

To expound on the resonance mechanism giving rise to the 
quench-independent frequency /„, it is crucial to note that, 
for an optimal quench J' = Jgp^, only the two singlet eigen¬ 
states IS' 1 . 2 ) vvith eigenvalues Ei ^2 are essentially involved in 
the dynamics, with (S'„|'!/i(0)) ~ 0 in Eq. (|6]) when n > 2. 
Namely, one numerically verifies that there is a small residual 
overlap 1 - I (^i IGS (ATi)) 1 2 - I (521GS (ATi)) 11 2 < 0.02 for any 
Ki> Kc and K 2 < when J' = J^p^. As one moves away 
from J' = Jgp^, other eigenstates rapidly come into play and 
significantly contribute to the time evolution in Eq. (|6]l. 

The dominance of two singlet eigenstates |5i_2) at the op¬ 
timal quench J' = J'^p^ suggests that the dynamics may be 
captured by an effective four-spin model. Consider 
and \\r)L/R for the impurity spins in the L/R parts, and 
I lt )L/fl and I l/r for the spins in the L/R bulks (see 
Fig. [ill- We represent |5i) ~ |0 “)l 0 |0“)_r, where |0^) = 
I t'lj') ± I iti)- Similarly, if we now make the ansatz that 

i>^2) — |1 )l <8 I — l)fl — |0 +)l 0 |0+)_R + I — 1)l 0 |1)_R, 

with |1) = |tlT) and | — 1) = |J,-IJ-), then the initial state 
|GS(A"i)) will be |GS(A:i)) ~ i|5i) - ^1^2). Thus, 
|T'(/)) ~ 5|5i) — e*^’^f'‘*^|52), which periodically brings 


N 

8 

12 

16 

20 

24 

28 

32 

36 

40 

fuN 

2.250 

2.184 

2.192 

2.180 

2.136 

2.100 

2.080 

2.01 

2.000 


TABLE I: The product f^N as a function of the length A: the table 
shows that fuN ~ 2.1 ±0.15 independently of the choice of the 
initial coupling Ki when K 2 < Kc- 


back |GS(Ari)). It is truly remarkable that the quench dy¬ 
namics of a complex quantum-many body system, when prop¬ 
erly tuned, can qualitatively be mimicked by four spins! In¬ 
tuitively, the entanglement makes the bulk spins collectively 
behave as two effective spins, forming a dynamically coordi¬ 
nated composite with the impurities. Putting this intuition on 
firm ground would be extremely interesting, and could open a 
new vista on quantum engineered quench dynamics. 


E. von Neumann entropy 

Since effectively one frequency governs the dynamics when 
J' = J'gpt, one may expect that the matrix elements of PL(t) 
oscillate with a few harmonics of fu- To verify this, we study 
the von Neumann entropy s{pL{t)) which depends on all lev¬ 
els of the entanglement spectrum. In Fig. |4ja), s{pL{t)) is 
given for two different quantum quenches. When increas¬ 
ing the difference K 2 — ATi, and thus, the energy released 
to the system through the quantum quench, deviations from 
a single frequency sinusoidal function become apparent, with 
the appearance of higher harmonics in the power spectrum 
S{f) = J^[s(t)]. This is shown in Fig. |4jb). A similar de¬ 
pendence on higher harmonics of the fundamental frequency 
fu is observed for the lower multiplets of the entanglement 
spectrum, i.e. for A„ with n > 2. 

Singlet fraction and spin correlations-- The emergence of a 
single frequency (and of its harmonics) is not a feature only 
of the dynamics of global quantities such as the entanglement 
spectrum or the von Neumann entropy. To probe for local 
quantities we trace out the bulks and compute the reduced den¬ 
sity matrix of the two impurities pii,,iR(f). This has the form 
of a Werner state due to the spin-rotational symmetry of the 
model, 

E \Tn{T^l (12) 

n=0,± 

where p{t) is the singlet fraction, \ip~) is the singlet and jT") 
are the triplets. The singlet fraction, which is nowadays ex¬ 
perimentally accessibl o^A33 ^ determines all local properties of 
the two-impurity composite, such as the two-point correlation 
= (1 ~ 4p)/3 and the two-impurity concurrence^ 
E{t) = max{0, 2p — 1}. In Fig.|4|c) the singlet fraction p{t) 
is plotted versus time for two different quantum quenches. As 
for the von Neumann entropy, the dynamics can be perfectly 
matched to the quench-independent frequency /„ for small 
quenches, and its higher harmonics (essentially the third one) 
for larger quenches. The Fourier transform P{f) = iF[p{t)] 
is plotted in Fig.Ud), showing the peaks for fu and 3/^. 
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FIG. 3: (a) The dominant frequency fpeak as a function of 1/J' in a semi-logarithmic plot for a chain of fixed length N = 20 when the 
impurity coupling J' varies. The deviation from linearity for J' < J'^p^ is due to the finite size effect as the entanglement length exceeds 
the system size, (b) The quench-independent frequency fu as a function of 1/ J'^pt on a semi-logarithmic plot. Each point corresponds to a 
different length N for which the optimal coupling J'^pt is found and then the quench-independent frequency /„ is determined through time 
evolution, (c) Data collapse for the dynamics of Ai(f) as a function of fut for three different lengths. 





Jt f 

FIG. 4: (a) The von Neumann entropy s{t) versus time for J' = 
J'opt = 0.3 in a chain of length N — 20 when K 2 = 0.19 and 
Ki = 1 (blue solid line) and Ki — 0.25 (red dashed line), (b) The 
Fourier transform S{f) of the von Neumann entropy versus / for the 
case of Ki — 1. The case for Ki = 0.25 (not shown in the figure) 
shows a similar peak at /„ with highly suppressed higher harmonics, 
(c) The singlet fraction p{t) versus time for J' = J'^pt = 0.3 in a 
chain of length N = 20 when K 2 — 0.19 and Ki = 1 (blue solid 
line) and Ki = 0.25 (red dashed line), (d) The Fourier transform 
P{f) of the singlet fraction versus / for the case of Ki — 1. 


IV. CONCLUSIONS 

In this article we have shown that a local quantum quench 
across the quantum critical point in a pertinently tuned two- 
impurity Kondo spin chain may lead to the emergence of long- 
lived single-frequency oscillations in the dynamics of the en¬ 
tanglement spectrum. The frequency thus revealed is indepen¬ 


dent from the energy released by the quench, and also shows 
scaling behavior with the system size, implying data collapse 
for the time evolution of levels in the entanglement spectrum. 
Important for possible experiments in the future, the quench- 
independent frequency leaves distinct hngerprints in local ob¬ 
servables, and it can be observed as the dominant frequency 
even when the system is not tuned to make the entanglement 
length extend over the full system. 

The fact that the single-frequency dynamics is found to 
be tied to the emergence of an optimal entanglement length 
may hint at new physics. By identifying the entanglement 
length with the dynamically generated screening length char¬ 
acteristic of Kondo systems^, our result points towards novel 
schemes for measuring it in the laboratory. Recent experi¬ 
mental progress in realizing two-impurity Kondo physics with 
tunable interactions^^ makes this a tangible challenge. 

We conjecture that a single frequency can be associated 
also to other massless systems. For this to happen there 
has to be a dynamically generated entanglement length 
which extends over the full system. In this work we used a 
spin-chain emulation of the two-impurity Kondo model as a 
paradigmatic example where this phenomenon occurs. Since 
the spin-chain emulator is a faithful realization of the spin 
sector of the two-impurity Kondo model, with a dynamics 
that effectively decouples from charge at low energies^, 
we expect that our results will be relevant also for the full 
two-impurity Kondo model with itinerant electrons. Thus, a 
spinful double-quantum-dot system, with the conhned spins 
interacting with the conduction electrons in their independent 
leads, could serve as an experimental setup where the emer¬ 
gence of a long-lived dynamics may be observed. 
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